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Weak radiative decay B —¥ X s "/ is known to be a loop-generated process. However, it does 
receive tree-level contributions from CKM-suppressed b — > uus'j transitions. In the present paper, 
we evaluate such contributions together with similar ones from the QCD penguin operators. For a 
low value of the photon energy cutoff Eo — mb/20 that has often been used in the literature, they 
can enhance the inclusive branching ratio by more than 10%. For Eo = 1.6 GeV or higher, the effect 
does not exceed 0.4%, which is due to phase-space suppression. Our perturbative results contain 
collinear logarithms that depend on the light quark masses m, (q — u,d,s). We have allowed 
mb/rn q to vary from 10 to 50, which corresponds to values of m q that are typical for the constituent 
quark masses. Such a rough method of estimation may be improved in the future with the help of 
fragmentation functions once the considered effects begin to matter in the overall error budget for 
B(B X.j). 

PACS numbers: 12.38.Bx, 13.20.He 



I. INTRODUCTION 

Weak radiative decay of the B meson is an invalu- 
able and well-established means for constraining physics 
beyond the Standard Model (SM). Its branching ratio 
has been measured to a few percent accuracy at the B- 
factories [l| • Theoretical calculations have acquired sim- 
ilar precision Q. The decay is generated dominantly by 
the b — > S7 transition that arises at one loop in the SM 
and its most popular extensions. However, it receives 
also tree-level contributions from the CKM-suppressed 
b — > uus^f transitions. In the present paper, we eval- 
uate such contributions together with similar ones that 
originate from the QCD penguin operators. 

It is well known that the inclusive B — > X s j decay rate 
can be approximated by its perturbative counterpart 



T(B -> X sl ) Ey>Eo ~ r(6 -> *J7K>E„. 



(1.1) 



where Xf stands for s, sg, sqq, . . . partonic states (with 
q = u,d, s only, as no charmed hadrons appear in X s 
by definition). Deviations from Eq. (jl.ip appear as cor- 
rections when Eq is large {Eq ~ hmb) but not too 
close to the endpoint (m& — 2E >• A ~ A QCD ). For 
Eq = 1.6 GeV, the corresponding non-perturbativc un- 
certainty amounts to around ±5% 0], while the known 
C(A 2 /ml c ) corrections 0,0] are smaller than that. 

The considered process is most conveniently analyzed 
in the framework of an effective low-energy theory which 
arises after integrating out the electroweak bosons and 
the top quark. The relevant effective weak interaction 
Lagrangian reads 



C 



weak 



4G F 



v; s v tb £ dPi + v: s v ub J2 CiiPi - p^) 



i=l 



(1.2) 

where Pi denote either dipole-type or four-quark oper- 
ators, and Ci stand for their Wilson coefficients. The 




FIG. 1: Tree-level Feynman diagrams for b — > sqqj at the LO. 
Black squares denote four-quark operator insertions. 



operators are given by 
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{sLl^l^l^T-bL^iqi^l^l^T-q), 

T ^m b (s L (T> iU b R )Fi ll/ , 

^m^SLa^T-b^G^. (1.3) 
^6 include all the active flavors q = 



Sums over q in P3 
U, d, s, c, b. 

Our goal in the present work is to evaluate O (a®) 
leading order (LO) contributions to the r.h.s. of Eq. (|1.1[) 
that originate from b — > qqsj with q = u,d, s. They can 
be generated either by the current-current operators P" 2 
or by the QCD penguin ones P3,...,6- The corresponding 
Feynman diagrams are shown in Fig. Q] 
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At first glance, it may seem surprising that the consid- 
ered LO effects have not been evaluated so far while the 
analysis of other contributions has already reached the 
O (a,) next-to-next-to-leading order (NNLO) level p|. 
Let us recall that the Wilson coefficients of all the op- 
erators in Eq. (| 1 . 3 j) acquire non-zero values already at 
the LO once the QCD logarithms have been resummed 
using renormalization group evolution from the elec- 
troweak scale (Xq ~ Myy,mt down to the low-energy 
scale /ib ~ mjj/2. However, despite being non- vanishing 
at the LO, Wilson coefficients of the QCD penguin 
operators remain rather small | C3 : ... : 6 (/^b)/C7 (Aib) | 2 < 
\Ci{pib)/C 7 {ji b )\ 2 ~ O.f, while the tiny CKM matrix el- 
ement ratio \V* s V u b/V t * s Vtb\ — 0.02 makes the current- 
current operators P" 2 even more suppressed. Moreover, 
when the lower photon energy cutoff Eq is at 1.6 GeV [2[ 
or higher the considered tree-level contributions to 
the branching ratio undergo severe phase-space suppres- 
sion, which justifies neglecting them at the leading and 
next-to-leading orders in a s (see App. E of Ref. @|). 

Given the current and expected future progress in the 
NNLO calculations Q, reliable uncertainty estimates in 
the SM prediction for B(B —> X s j) can no longer be 
made without evaluating the diagrams in Fig. [T] and 
checking what the actual size of their contribution to the 
r.h.s of Eq. (jl.ip is. This fact serves as the main mo- 
tivation for our present work. Since the corrections are 
expected (and found) to be quite small for Eq = 1.6 GeV, 
rough estimates of their size are sufficient. Actually, 
nothing more is available within perturbation theory 
alone because collinear logarithms \n.(mb/m q ) involving 
light quark masses m q (q — u, d, s) remain in the final 
expressions. While using the so-called current masses for 
the light quarks is not adequate in such a case, the above- 
mentioned rough estimates can be obtained by assuming 
that m q are of the same order as masses of pions and 
kaons or, equivalently, as the constituent quark masses. 



We shall do it by varying mb/m q from 10 to 50, which 
covers the necessary range. A refined approach would 
require taking non-perturbative fragmentation into ac- 
count, as it has been done in Refs. [1, Q for contribu- 
tions that are proportional to \C%\ or in Ref. [13 for the 
b — > uudj background. Such an analysis is beyond the 
scope of the present paper. 

The article is organized as follows. Our final perturba- 
tive results together with a discussion of their numerical 
impact on the total decay rate are presented in Sec. [TXJ 
The next two sections contain brief descriptions of two 
alternative methods that we have used for integration 
over the four-body phase space. In Sec. IIII1 partly 
massive phase-space integration in D = 4 dimensions is 
outlined. In Sec. IIV1 a calculation involving dimensional 
regularization and splitting functions is described. We 
conclude in Sec. [Vj App. A contains a collection of 
several intermediate results that may be useful for 
other calculations of electromagnetic bremsstrahlung 
corrections to decays mediated by four-fermion opera- 
tors. App. B is devoted to summarizing the necessary 
properties of the splitting functions. 



II. RESULTS 

Our final result can be expressed in terms of three 
functions T&(<5) (k — 1,2,3) that depend on the pho- 
ton energy cut Eq = ^(1 — 8) and on logarithms of the 
quark masses. Each of the functions gets multiplied by a 
quadratic polynomial in the Wilson coefficients Cj that 
are evaluated at the scale We assume that all the Cj 
are real, as it is the case in the SM. The LO tree-level 
contribution to T[b — > Xfj] arising from the four-quark 
operators P" 2 an d ^\...,6 reads 



AT[b -> Xfj]^ = 
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where Ai = —Ci Xy?Y? b , i = 1, 2, and 
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Ci = -0.8144 


Cz = -0.0125 


C 5 = 0.0012 


CV = -0.3688 


C 2 = 1.0611 


C 4 = -0.1224 


C 6 = 0.0026 


C 8 = -0.1710 



TABLE I: The LO Wilson coefficients C % at fi b = 2.5 GeV. The 
matching scale fio has been set to 160 GeV in their evaluation. 
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P (S) = 5+-5 4 + \n(l-S), 
6 

u(S) = ^S 2 -2S 3 + S\ 



(2.5) 



The function Ts(d) originates from cross-terms in b — > 
sss^f where the s-quark lines are interchanged in one of 
the interfered diagrams (see Fig. [2]). All the other con- 
tributions from the penguin operators f3,....6 alone are 
described by Ti(8). Finally, T 2 (S) comes from P™ 2 and 
their interference with P^.,.,6- 

We have retained the light quark masses m q in the 
collincar logarithms only, i.e., all the power-like correc- 
tions proportional to m 2 / mf have been neglected in the 
above expressions. Such an approximation breaks down 
at some point, which manifests itself in non-physical 
negative values of Ti^ 2 (S) when ln(rn^ S j 'to 2 ,) is not big 
enough. 

Determining the size of the calculated correction is 
now straightforward. Numerical values of the LO Wil- 
son coefficients Q = C\ ■ (pb) are summarized in Tab.|U 
For the CKM element ratio we use XhY? b = -0.0079 + 

0.018 i 11]. As far as the light quark masses m Ut d, s are 
concerned, we set all of them equal in the numerical ex- 
amples to be discussed below. 

In Tab. UH we present the calculated correction (|2.1[) 
as a fraction of the leading contribution to the decay 
rate = G 2 F a cm m 5 b C% \V t * s V tb \ 2 /(32tt 4 ) for two val- 
ues of Eg and two values of ^+ Strong dependence on 
the collinear logarithms is clearly visible. On the other 
hand, the non-logarithmic terms turn out to be relevant 




FIG. 2: Sample contribution to the b — > sssj cross-terms that 
give rise to the function Ts(8). It is represented as a cut propa- 
gator diagram, with the cut denoted by the vertical dashed line. 



TABLE II: The correction from Eq. (|2.ip as a fraction of the 
leading contribution r' ' for several values of Eq and mb/m q . 



in the considered range of ■j-j^. For a low value of the 
photon energy cutoff Eq ~ m(,/20 that has often been 
used in the literature, the correction can enhance the 
inclusive rate by more than 10%. On the other hand, 
for Eo = 1.6 GeV ~ ^ or higher, the effect does not 
exceed 0.4%, which is obviously due to phase-space sup- 
pression that becomes efficient when we approach the 
high-energy endpoint Eq ~ Sit., i.e. when S becomes 
small. In this limit, our correction in Eq. (12.11) behaves 
like O (f5 2 In S) . In all the four cases shown in Tab. [Til the 
contribution from T\(S) is the dominant one, while T^IS) 
(T 2 (<5)) give about 5 (50) times smaller effects. Thus, the 
CKM-suppressed contributions that come with T 2 (8) are 
minuscule indeed, and no precise knowledge of V y* y-" b is 
necessary here. On the other hand, our LO results ex- 
hibit significant dependence on the renormalization scale 
jj,b that comes from C^,. ..^(Mb)- It could be stabilized 
only after including 0(a s ) contributions to b — > qqsj, in 
particular the ones generated by Pi,2- 



III. CALCULATION INVOLVING PARTLY 
MASSIVE PHASE-SPACE INTEGRALS 

Let us now briefly describe the calculation. First, we 
consider the diagrams in Fig. Q] with an operator 



Pa = (s L j fi b L )(q 1 YQ2) 



(3.1) 



where no sum over flavors is present (contrary to P3 in 
Eq. (II. 3|) ). and the electric charges {Q s , Qb, Qi, Q2} are 
retained arbitrary. The invariant matrix element M. is 
calculated in the Feynman gauge, so collinear divergences 
are allowed to occur in interferences between different di- 
agrams (rather than in self- interference terms alone) [l2j . 
The Dirac algebra is performed in D = 4 — 2e dimensions 
without neglecting the light quark masses. At this point, 
we proceed in two different ways. One of them is to in- 
tegrate the spin- averaged |WJ| 2 over the partly massive 
four-body phase space in D = 4 dimensions. The other 
way is to initially neglect the light quark masses, inte- 
grate \ftA\ 2 in D — 4 — 2e dimensions over the massless 
phase space, and convert the collinear 1/e poles to loga- 
rithms of masses only afterwards (see Sec. HV)) . 

In the D = 4 case, we impose the photon energy cut (in 
the 6-quark rest frame J-) after expressing P 7 in terms of 



invariants, namely 2mi>E 1 



S123, where S123 is the 



invariant mass squared of the q\q 2 s system. A boost to 
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the rest frame T of this system is performed along the — k 
direction, where k is the photon three-momentum in J-. 
The qi&s system has energy equal to v /si23 in J 7 , while 
the three-momenta of its constituents define a plane. The 
direction of k with respect to this plane is parametrized 
by two polar angles 9 and if. The remaining three phase- 
space variables are the quark energy fractions Xj — ^7== 
(i = 1,2) and S123. 

We integrate first over 9 and <p, and obtain results 

2 

containing large logarithms In — | . Integration of the in- 



terference terms 



involves angular ordering (see 



Ref. [13]), which applies also to the case \r • ,J 



m q ^ m b- 

At this stage, all the collinear logarithms have been al- 
ready identified, which allows us to neglect masses of the 
outgoing particles in the remaining terms. Next, integra- 
tions over the quark energy fractions and S123 are per- 
formed. Finally, charge conservation Qb = Qi — Q2 + Qs 
is imposed, but the charges on the r.h.s. are still re- 
tained arbitrary. Several intermediate results with such 
arbitrary charges are collected in App. A. 

Once the final result for P3 with arbitrary charges is 
at hand, obtaining the corresponding one for P3 (T\(S)) 
is just a matter of substituting the actual values of the 
charges and summing over flavors. Extending the calcu- 
lation to ^4,5,6 and taking into account the b — > sssj 
cross-terms (T 3 (6)) requires modifying the Dirac alge- 
bra and color factors but no essential difference in the 
phase-space integration is encountered. As far as 72(5) 
is concerned, it originates from the operators P" 2 , in- 
cluding their interference with Pa,. ..,6- In this case, it is 
sufficient to express P" 2 in D = 4 dimensions as linear 
combinations of the (s6)(uu)-parts of P3,...,6, namely 



pit 4 rtu 

-"1 — 27*3 



1 T)U 1 1 nil 1 JDU 

5/4 i- 2f-r 5 3g-r 6 



pit 1 jdu 2 r>u 1 1 pu 

r 2 — 9^3 3^4 T 36^5 



1 pit 
6^6 



(3.2) 



The above expressions are easily derived using Fierz iden- 
tities for the Dirac and Gell-Mann matrices. 



IV. CALCULATION WITH THE HELP OF 
DIMENSIONAL REGULARIZATION 

A calculation with the help of dimensional regulariza- 
tion is technically simpler but involves a few subtleties. 
To start with, all the particles in the final state are as- 
sumed to be massless, and the phase-space integration is 
performed in D — 4 — 2e dimensions [14|. The results 
are collected in App. A. Given that the collinear diver- 
gences appear as 1/e poles, one should make sure that 
no ambiguities arise from Dirac traces with odd numbers 
of 75's [15]. Fortunately, such traces being purely imagi- 
nary give no contribution to the decay rate in the case of 
P3,... 6 alone. As far as P" 2 are concerned, we do not need 
to consider them at this level. Once the collinear diver- 
gences are re-expressed in terms of logarithms of masses 



(see below), we can pass to D = 4 and use the identities 
(& 

For definiteness, let us consider the operator P 3 from 
Eq. (|3.1j) again, but this time with the s-quark denoted 
by 53 to keep the notation symmetric. Before integrating 
over the photon energy (but after integration over all the 
other phase-space variables) , the differential decay width 
for b — > gig2937 reads 



rfT _ cfT, fir shift 
dx dx dx 



(4.1) 



where x = 2Ej/mb- The first term on the r.h.s. above is 
the dimensionally regulated expression, while the second 
one converts the dimensional regulators to logarithms of 
masses. Its explicit form is given below. 

We shall need a D-dimensional expression for the to- 
tal width of the three-body decay b —> qiq 2 q3- De- 
noting momenta of the final-state massless particles by 
Pi, i = 1,2,3, and parametrizing the phase space by 
sij = 2piPj/ml, one can write 

r 3 -body = 7T— / ds i2 ds i3 ds 23 \M'\ 2 M 3 , (4.2) 

where Ai' is the corresponding invariant matrix element, 
J1 2 = ^ 2 e 7E /47r, and 

2^T(l-e)T(l-e)( Sl2Sl3S23 y 1 '*> 

describes the phase-space measure [lj| • The second term 
in Eq. (14. ip can now be written as 



dx 



2m;, 



d Sl2 ds 13 ds 23 \M'\ 2 M 3 ^ 



- 1 + 2 In 



1 - S23 
(1 - s 23 )ii 



2ttx 

e(i-s 23 -z) 

+ (cyclic) I . (4.4) 



Properties of the splitting functions that have been nec- 
essary to derive the above formula are summarized in 
App. B. One should remember that all the collinear 1/e 
poles cancel in Eq. (14. ip only after the charge conserva- 
tion Qb = Qi — Q 2 + Q s has been imposed. 

The structure of Eq. I|4.4p remains the same irrespec- 
tive of what interaction generates the 6-quark decay. 

Thus, it is applicable as it stands to the operators P3 6- 

Next, as already mentioned, Eq. (|3.2I) is used in D = 4 
dimensions to take P" 2 into account. This way the final 
expression in Eq. (I2.lt has been obtained once again, in 
perfect agreement with the results of Sec. IIIII 



V. CONCLUSIONS 

In the present paper, we have evaluated the LO con- 
tributions to the partonic decay width T(b — > Xfj) that 
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originate from the four-quark operators P" 2 and p3....,6- 
They can be sizeable (above 1 0%) for low photon energy 
cutoffs Eq but become very small (below 0.4%) in the 
phenomenologically interesting domain Eq > 1.6 GeV, 
i.e., for 6 = 1 — 2E^/m\ } < 0.32. For small 8, they be- 
have like 0(S 2 \n8), which determines their phase-space 
suppression near the endpoint. 

The presence of collinear logarithms \n(rn 2 S/m 2 ) in- 
volving light quark masses m q (q — u, d, s) implies 
that our perturbative results (with mb/m q varied from 
10 to 50) may serve only as rough estimates of the 
corresponding contributions to the inclusive hadronic 
B — > X s "f decay rate. Having such rough estimates 
at hand is both advantageous and sufficient at present, 
given that the overall non-perturbative uncertainty re- 
mains at the ±5% level @. However, once our control 
over non-perturbative corrections improves in the future, 
the current contributions will need to be supplemented 
with hadronic fragmentation effects, along the lines of 
Refs. [li]. 
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APPENDIX A: INTERMEDIATE RESULTS 

Here, we present several intermediate results that 
might be useful for studying hard photon emission in 
other processes mediated by four-fermion operators, like 
muon decays or semileptonic heavy quark decays. Al- 
though radiative corrections to such processes have been 
calculated long ago, none of the published results that 
we are aware of leaves the photon energy as the only 
phase-space variable that has not been integrated over 
(or, equivalently, integrated with an arbitrary cutoff). To 
make our results applicable to such cases, it is enough to 
present them for arbitrary electric charges of the final- 
state fermions. 

As in Sec. |Hll we replace P 3 ([O]) by P 3 flHH]) and 
assume that q\ ^ s, which means that no cross-terms 
(T 3 {8)) arise. The contribution of P 3 to the decay rate 
is as in the Cf Ti(5) term in Eq. (|2~Tj) but with T X (S) 



replaced by 

fx (6) = {Ql + Q 2 2 ) Fn(S) + Q 2 S F SS (S) 
+ QiQ2(-\ P {S)-\u{5)-\5 2 



+ (Q\ -Q-i)Q, ( \ P {5) + ^uj{5) j . (A.i) 
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F ss (6) = -p(S) In- 



—s - —S 3 + —8 A 
6 12 144 



17 
7T 



ln(l - 5) - Li 2 (<5) 



(A.2) 



In Eq. (|A.1|) . charge conservation Qb = Qi — Qi + Q s 
has been already imposed. In effect, collinear logarithms 
remain only in the terms that come with Q 2 , i = 1,2, s. 
For simplicity, all the light quark masses have been set 
equal and denoted by m q . However, it is easy to relax 
this assumption and identify them by the corresponding 
charges despite the fact that charge conservation has al- 
ready been used (see App. B). 

Now, let us consider the case when qi = s in Eq. (13.11) . 
Then, apart from Ti(5) (jA.ll) . we get an additional con- 
tribution from the cross-terms 

T 3 (6) = Q 2 2 S 22 (5) + Q 2 S S SS (S) 

+ Q2Q S (™p(*) - ^oj(S) - M , (A.3) 



where 

_ Us 3 + —S 4 + — ln(l - 5) - -LUS) , 
12 144 18 V ' 3 W ' 

s ss (s) = - 2 - P (s) +^ s -r9 s2 -r^ 3 

3 m q 18 12 18 

70 A1. 9 

+ 216 54+ 18 ln(1 - 5) -3 Li2W - (A ' 4) 

The above results are valid only for the particular Dirac 
structure of the operator P3. To generalize them to all 
the four-fermion operators with chirality-conserving cur- 
rents, it is sufficient to consider 

h = (sLl^ 2 lM(q ll ^r 2 r 3 q2)- (A.5) 

Its interference with P3 gives 
20T 1 (5) + (Ql - Ql) n{8) + {Q 1 + Q 2 ) Q s r 2 {5), (A.6) 
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and 32T3((5) for the cross-terms. Its self-interference 
gives 

136Ti(<5) + 10 [{Q\ - Ql) n(S) + (Q a + Q 2 ) Q s r 2 (6)] , 

(A.7) 

and 256 T^J) for the cross-terms. The functions ri.2(<5) 
are given by 

m 2 A 31 
n(S) = 2cj(<S)ln— - 2p{5) - IIS 2 + 16S 3 - — S 4 , 



r 2 {8) = Ap{8) - -u{S) - 3<5 2 . 



(A. 



The corresponding results for massless hnal-state 
fermions in the dimensional regularization have been ob- 
tained with the help of Eq. (A. 5) of Ref. [16j where a com- 
pact expression for the Z)-dimensional four-body phase- 
space measure is given. Only the functions that come 
proportional to squared charges differ from their D = 4 
counterparts. For P3 alone, they read 



for the P5 self-interference (Eq. (|A.7[) and below) read 

(Ql-Qt)n(5) 



2(Q 2 -Q 2 ) 



C +^ 



S 2 

iop- t 



Ql 



■18p+ ~L) 



256 T a 



256 T. 



(<0 



64 
~9 



4Q 2 [9p+^uj 



32 



19 
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Q;(21p-4w)-yQ2(7p+5w). 

(A.ll) 



Conversion of the collinear regulators in Eq. (|4.ip is most 
conveniently performed before substitution of charges, at 
the level of the functions from Eqs. (|A.9I) — (| A. 1 1|) . 

Our final result in Eq. (|2.1[) has been obtained by form- 
ing appropriate linear combinations of the functions ap- 
pearing in Eqs. (|A.1[) . (|A.3|) and (|A.8[) . according to the 
values of charges, color factors and sums over flavors. 



APPENDIX B: SPLITTING FUNCTIONS 



F$(5) = + I w (*))c B (fl + ^+|a*-|«» 

+ ^ + ^ ln(l -S) + ln 2 (l - S) 2U 2 (8) , 

F$(S) = p(6)C t (S)+ 1 -§S+^S 2 + ^ + ^ 

+ ^ ln(l -S) + ln 2 (l -6)- 2Li 2 (<5) , 



44 r3 83 r4 101. „ „ 1, a/, c\ 
- §Li 2 (5), 

+ ^\n(l-5) + hn 2 (l-5)-^U 2 (5), (A.9) 

where C e (6) = l/e-21n[5(l - 5)}. 

In the case of P3-P5 interference (Eq. (IA.6I) and below) , 
the following replacements need to be made: 



(Qi - Q2) n -2(Qi-Qi) \(C e + 8)u + 7p-- 
-UQlp-AQliGp + u), 
32 T 3 -> 32T 3 (e) -^Q2 (3/9 _ w) _8 g 2 (/g + w); 

(A.10) 

where (5) is given in terms of (6) and sil* (S) , in 
analogy to Eq. (|A.3[) . The corresponding replacements 



Our conversion formula (|4.4[) involves a difference be- 
tween splitting functions derived in the dimensional reg- 
ularization and in the regularization with masses. Their 
derivation along the lines of Refs. [l?], EH * s briefly de- 
scribed in the following. 

Let us consider an amplitude M.(q, k; . . .) of a process 
where an external massive fermion radiates a photon with 
momentum k (k 2 = 0). After the radiation, the fermion 
is on shell (q 2 = to 2 ). We shall assume that the calcula- 
tion is performed in the light-cone axial gauge n ■ A = 0, 
where n is a lightlike vector (n 2 = 0), and the sum over 
photon polarizations gives 



y- a* a 



kn 



(B.l) 



In such a gauge, interference terms between different di- 
agrams are free of collinear singularities fl2l ]. 

In the course of the splitting function derivation, it is 
convenient to introduce the Sudakov parametrization in 
terms of 

qk 



P 



k 



(q + k)n 



k X = zp - q + z(l - z) ( q + k > q n , (B.2) 

qn 

where z = (qn) / [(q + k)n] <G [to 2 /(to 2 + 2qk) , l] . It is 
easy to verify that p 2 = m 2 , and that the spacelike vector 
k± is orthogonal to both n and p. Inverting the relations 
Eg) , one finds 



q = zp — k_y 



k'j + (1 - z 2 )m 2 
z (2pn) 



n, 



k 2 - (1 - z) 2 to 2 , , 
k = (l-z)p + k ± + ^ Uo \ n, B.3) 

(1 - z) (2pn) 
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The integration is performed from = up to — 
E 2 , where E is chosen to remain in the \ow-k± region, to 
preserve the factorization formula. 

The splitting functions integrated over k^_ read 



FIG. 3: A schematic picture of the factorization formula (|B.5 



Mz) = - 

7T 



l+Z 



1 , E 

z \ 2e jj. 



1 



where k^ stands for —k^k± ll when expressed in a frame- 
independent way. In such a parametrization, propagator 
denominators that are responsible for collinear singular- 
ities appear as 

1 1 z(l-z) 

a 2 2qk 



ki 



+ (1 — z) 2 m 2 



(B.4) 



The transverse momentum k± parametrizes how far off- 
shell the radiating fermion is. In the massless case (m = 
and e 7^ 0), the collinear limit is defined by k± — > 0, which 
determines the phase-space region where the 1/e singu- 
larity arises. In the case of a massive fermion (m ^ and 
e — 0), the quasi-collinear limit has to be considered (l9j . 
In this limit, the collinear region is defined by taking si- 
multaneously kj_ — > and m — > 0, but keeping the ratio 
m 2 /k 2 , fixed. Both limits lead to the factorization for- 
mula fnl US] illustrated in Fig. |3] 



\M{q,k-...){ 



^P(z)\M(p;...)\ 2 . 



(B.5) 



where A4(p; . . .) is the amplitude of the process with- 
out radiation where p — > q/z in the collinear or quasi- 
collinear limits, and Qj is the fermion charge. The split- 
ting function P[z) in the collinear limit and in D = 4— 2e 
dimensions reads (see Eqs. (53) and (54) of Ref. fl7| ) 



P t {z) = Sir a e 



1 



1 



(B.6) 



It becomes the Altarelli-Parisi splitting function [2l| for 
the gluon emission off quark when 8Tra e is replaced by 
Cp. Its extension to the massive quark case in D = 4 
dimensions is 



P m (z) = 8i:a e 



f 



1 



qk 



(B.7) 



For definiteness, let us assume that the number of final- 
state particles is as in Fig [TJ In the collinear region, it 
is possible to disentangle the four-body phase space d$4 
into a convolution of the three-body phase space of the 
non-radiative process, and the pha se space corresponding 
to the radiation process alone [20L |22| , 

d$4 = d$ 3 ® tf$ . (B.8) 
One proceeds with integration over the low-fc_L region 
using the following phase-space measure [17J 

d d k 

d T S + (k 2 ) = 



(A) 



(27r) d " 
1 1 
16tt 2 T(l - e) 



dkj 



r,2 



dz 



e(z(i-z)). 

(B.9) 



fm(z) = 

7T 



1 + z 2 

I- z 



In 



E 



(1 — z)m 1 — z 



, (B.10) 



while their difference is 

Af(z) = f m (z)-f e (z) 
a e 1 + z 2 r 



2?r 1 - z 



l-l-21n (1 ~ z)m 
e n 



(B.ll) 



The dependence on E cancels in Eq. (|B.11|) because both 
splitting functions in Eq. (IB. 10[) have been consistently 
derived in the corresponding regularizations, and they 
contain the same high-fcj^ finite terms. The formula 
(jB.lip could have also been obtained with the splitting 
functions from Ref. [JJ . 

The mass-regulated and dimensionally regulated ra- 
diative decay widths with E~ t > Eq satisfy the following 
relation 



z)\M(p;...t 



Q[(l -z)p° j -E ] 



(B.12) 



where the sum goes over all the radiating fermions. 
The three-body final-state phase-space measure rf <&3 in- 
tegrated over the "blind" angular variables reads [16[ 



d® 3 = /T £ M 3 dsi2 ds 13 ds 2 



(B.13) 



where M 3 has been given in Eq. (14.3[) . 

Eq. (|B.12|) has been derived in the light-cone axial 
gauge but it is actually gauge-independent once charge 
conservation has been imposed. Thanks to this fact, we 
could have used it in our Feynman-gauge calculation. 
Actually, the conversion formula (|4.4[) is obtained from 
Eq. (|B.12|) just by differentiation with respect to Bo- 
Even after imposing charge conservation, the splitting 
functions derived in the same regularization may differ 
by finite terms which depend on the chosen gauge and 
also on the high-fc^ integration limit. Only the differ- 
ence (|B.11[) of the two splitting functions is gauge- and 
convention- independent . 
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